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High precision, absolute total column ozone measurements 
from the Pandora spectrometer system: Comparisons with data 
from a Brewer double monochromator and Aura OMI 

12 23 24 25 

Maria Tzortziou, Jay R. Herman, ’ Alexander Cede, ’ and Nader Abuhassan ’ 

Received 22 March 2012; revised 15 June 2012; accepted 3 July 2012; published 22 August 2012. 

[i] We present new, high precision, high temporal resolution measurements of total 
column ozone (TCO) amounts derived from ground-based direct-sun irradiance 
measurements using our recently deployed Pandora single-grating spectrometers. 

Pandora’s small size and portability allow deployment at multiple sites within an urban 
air-shed and development of a ground-based monitoring network for studying small-scale 
atmospheric dynamics, spatial heterogeneities in trace gas distribution, local pollution 
conditions, photochemical processes and interdependencies of ozone and its major 
precursors. Results are shown for four mid- to high-latitude sites where different Pandora 
instruments were used. Comparisons with a well calibrated double-grating Brewer 
spectrometer over a period of more than a year in Greenbelt MD showed excellent 
agreement and a small bias of approximately 2 DU ( or, 0.6%). This was constant with slant 
column ozone amount over the full range of observed solar zenith angles (15-80°), 
indicating adequate Pandora stray light correction. A small (1-2%) seasonal difference was 
found, consistent with sensitivity studies showing that the Pandora spectral fitting TCO 
retrieval has a temperature dependence of 1% per 3°K, with an underestimation in 
temperature (e.g., during summer) resulting in an underestimation of TCO. Pandora agreed 
well with Aura-OMI (Ozone Measuring Instrument) satellite data, with average residuals 
of <1% at the different sites when the OMI view was within 50 km from the Pandora 
location and OMI-measured cloud fraction was <0.2. The frequent and continuous 
measurements by Pandora revealed significant short-term (hourly) temporal changes in 
TCO, not possible to capture by sun-synchronous satellites, such as OMI, alone. 

Citation: Tzortziou, M., J. R. Hennan, A. Cede, and N. Abuhassan (2012), High precision, absolute total column ozone 
measurements from the Pandora spectrometer system: Comparisons with data from a Brewer double monochromator and Aura 
OMI, J. Geophvs. Res., 117, D16303, doi:10.1029/2012JD017814. 

1. Introduction 

[ 2 ] The importance of ozone (O3) in Earth’s stratospheric 
and tropospheric chemistry, and its influence on radiative 
forcing and climate change, surface and underwater UV-B 
levels, local and regional environmental degradation, human 
health and vegetation, are well known, and have been 
discussed extensively in the literature [Levy, 1971; Fishman 
et al., 1979; Crutzen, 1995; Carpenter et al., 1997; Heck 
et al, 1982; Reich and Amundson, 1985; Diffey, 1991; 
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Tevini, 1993; Chameides et al., 1994; Environmental 
Protection Agency, 1998; World Health Organization, 2000; 
World Meteorological Organization, 2006; Intergovernmental 
Panel on Climate Change, 2007; Zerefos et al, 2009; Hennan, 
2010]. Meteorology and photochemical mechanisms involv- 
ing various species of natural and anthropogenic sources (e.g., 
NO x , CO, CH4, VOCs, CFCs) have been shown to influence 
the formation, dispersion, destruction, and transport of O 3 in 
the atmosphere, often resulting in significant heterogeneities 
in O3 distribution [ Fishman et al., 1979; Comrie and Yamal, 
1992; Solomon, 1999; Tu et al, 2007; Duehas et al, 2002; 
Lehman et al, 2004; Crutzen, 1995; Crutzen et al, 1999; 
Chameides, 1978; David and Nair, 2011]. Nitrogen dioxide 
(NO 2 ) plays a critical role in these processes, as a major O 3 
precursor in the lower troposphere and through catalytic 
destruction at higher altitudes [ Chameides , 1978; Logan et al, 
1981; Brasseur et al, 1998]. Simultaneous measurements of 
O3 and NO2 are, therefore, essential for improved under- 
standing and predictive modeling of ozone photochemistry, 
distribution, and dynamics. Toward this objective, new 
ground-based spectrometer systems, Pandora, were recently 
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deployed at multiple sites in the U.S. (12), Europe (2) and 
Korea (2) for high temporal resolution, high precision, simul- 
taneous measurements of various trace gases including 0 ; , 
N0 2 , S0 2 , BrO, water vapor and formaldehyde. Total column 
N0 2 retrievals from Pandora have been discussed in Herman 
et ah [2009] and shown to be of high precision (0.01 DU) 
and good accuracy (±0.05 DU). The purpose of this paper is to 
validate the capability of the Pandora spectrometer system to 
measure total column ozone (TCO) amounts compared to the 
widely accepted Brewer double monochromator throughout 
the day, and to the AURA-OMI satellite instrument during its 
overpass time. 

[3] Ground-based measurements of TCO have been made 
during the past decades using a wide variety of instruments, 
including filter instruments and spectrometers. Dobson and 
Brewer spectrometers [Dobson, 1957; Kerr et al., 1981; 
1988] have been the most widely used, as they are the 
standard instruments for the Global Atmosphere Watch 
(GAW) program of the World Meteorological Organization, 
and remain one of the main sources of information about 
long-tenn ozone changes, trend analyses, and identification 
of possible natural or anthropogenic influences on atmo- 
spheric composition [Loyola et al., 2009; Fioletov et al., 
2011]. In addition to TCO retrievals, Brewer measurements 
have been applied to retrievals of UV spectral measurements 
[e.g., Lakkala et al., 2008], aerosol optical depth [e.g., De 
Bock et al., 2010], total column S0 2 and N0 2 amounts 
[Kerr et al., 1988; Fioletov et al., 1998; Cede et al., 2006b] 
and O3 vertical distribution [Tzortziou et ah, 2008]. Dobson 
(starting in mid- 1920s) and Brewer (developed in the late 
1970s) instruments have a long history of ozone measure- 
ments that spans the data record for satellite retrievals 
(1979 to present) from TOMS, SBUV, SBUV-2, GOME, 
SCIAMACHY and OMI (Total Ozone Mapping Spec- 
trometer, Solar Backscatter Ultraviolet, Global Ozone 
Monitoring Experiment, SCanning Imaging Absorption 
spectroMeter for Atmospheric CHartographY, and Ozone 
Monitoring Instrument) [Bhartia et al, 1996, Burrows et al., 
1999, Bovensmann etal., 1999; Veefkind and deHaan, 2002; 
McPeters et ah, 2008]. TCO retrievals for all instruments 
discussed here are based on the spectral fitting technique [see 
e.g., Noxon, 1975; Platt et ah, 1979], where the laboratory- 
measured absorption features of different species are fitted in 
the measured signals. Older instruments, such as Brewer, 
Dobson, TOMS, or SBUV use <10 channels (wavelengths), 
while newer instruments such as Pandora, SCIAMACHY 
and OMI use a much larger number of wavelengths in their 
retrievals. The Pandora TCO retrievals discussed in this 
paper use a variance weighted spectral fitting algorithm on 
the Pandora ~0.6 nm spectral resolution direct-sun radiance 
measurements in the spectral range 310-330 nm, after cor- 
rection for stray light effects. 

[ 4 ] Improved understanding of small-scale atmospheric 
dynamics and spatial heterogeneities in trace gas distribution 
through ground-based observations require measurements at 
high spatial and temporal resolution. The small size and 
portability of the Pandora spectrometer system allow 
deployment at multiple sites within an urban air-shed, and 
development of a ground-based monitoring network for 
studying trace gas dynamics, local pollution conditions, 
photochemical processes and interdependencies of ozone 
and its major precursors. Such information is also critical for 


validation of satellite trace gas retrievals under a wide vari- 
ety of conditions [Herman et ah, 2009]. 

[5] In this paper we present high temporal resolution 
measurements of absolute total column ozone amounts 
derived from the Pandora spectrometer. For the first time, 
ozone retrievals are routinely performed using a high- 
resolution extraterrestrial spectrum derived from indepen- 
dent measurements as the reference spectrum. The spectrum 
used is the Kurucz [2005] spectrum that is radiometrically 
corrected to the Atlas-3 SUSIM spectrum (M. E. Van Hoosier, 
The ATLAS-3 solar spectrum, 1996, available at ftp://susim. 
nrl.navy.mil/pub/atlas3) as described in Bernhard et ah 
[2004]. Convolution of the Kurucz spectrum to the Atlas-3 
SUSIM resolution enables the removal of a small radio- 
metric variation in the Kurucz spectrum that could mas- 
querade as an atmospheric absorber if uncorrected. The 
result is a radiometrically corrected Kurucz spectrum with 
very high wavelength accuracy (picometer) and an estimated 
radiometric offset accuracy of 3% [Thuillier et ah, 2003]. 
The corrected high-resolution spectrum was then convoluted 
with the Pandora slit function to be used as a top-of-the 
atmosphere (TOA) reference spectrum. This method is very 
different from the usual technique of using a TOA reference 
spectrum based on the instrument’s own measurements from 
a Langley extrapolation to zero air mass [e.g., Herman et ah, 

2009] . Using the independent spectrum technique gives an 
absolute measurement of TCO, since the reference spectrum 
is free of residual atmospheric absorption. The method 
requires an accurate knowledge of the instrument’s slit 
function as determined in the laboratory, but is insensitive to 
a radiometric calibration offset that only changes smoothly 
with wavelength. For example, a change that can be fit with 
a 4th order polynomial. 

[6] Results are presented for four mid- (i.e., 39°N) to high- 
(i.e., 65°N) latitude sites where different Pandora instruments 
(one CMOS system #9, and two CCD systems #2 and #21) 
were used: GSFC (Goddard Space Flight Center), in Mary- 
land U.S., where the longest Brewer and Pandora collocated 
data records have been available (June 2009 to September 

2010) ; Cabauw in the Netherlands, where measurements were 
obtained as part of the CINDI (Cabauw Intercomparison 
Campaign of Nitrogen Dioxide measuring Instruments) cam- 
paign during June-July 2009 [Roscoe etal., 2010]; Helsinki in 
Finland (September-October 201 1); and Fairbanks in Alaska, 
where observations were performed as part of a large field 
campaign focusing on trace gas measurements at high latitudes 
(March-April 2011). Comparisons between Pandora, our 
ozone network double-grating Brewer spectrometer (#171), 
and Aura-OMI satellite retrievals are discussed. 

2. Methods 

2.1. Pandora Technical Characteristics 

[ 7 ] The Pandora spectrometer system consists of a head 
sensor with foreoptics, mounted on a computer controlled 
sun-tracker and sky-scanner (^0.01° pointing precision), 
and connected to an array spectrometer by means of a 400 
micron single strand multimode optical fiber (Figure 1). A 
drawing of the head sensor is shown in Figure 2. To achieve 
wavelength and radiometric calibration stability, the spec- 
trometer is temperature stabilized inside an insulated enclo- 
sure using an actively coupled thermo-electric cooler and 
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Figure 1. Sky-radiance plus direct-sun Pandora CCD spec- 
trometer system (left) and direct-sun only Pandora CMOS 
spectrometer (right) in Cabauw. 

heater. A detailed description is given in Herman et al. 
[2009], 

[8] There are two versions of the Pandora spectrometer 
system, the Sun-only CMOS detector Pandora and the CCD 
detector Sun-and-Sky Pandora (Table 1). The CMOS Pan- 
dora spectrometer system operates only in direct-sun 


Table 1. Sun-Only CMOS and Sun-and-Sky CCD Detector 
Pandora Characteristics 



Sun-Only CMOS 

Sun-and-Sky CCD 

Wavelength Interval 

280 to 500 nm 

280 to 525 nm 

Spectral resolution 

0.6 nm 

0.6 nm 

Oversampling 

2.3 pixels 

4.2 pixels 

Detector 

CMOS 1024 pixels 

2048 pixels 
backthinned CCD 

Field of view 

1.6° FWHM 

1.6°FWHM 

Pointing Precision 

0.01° 

0.01° 

Neutral Density Filter 

None 

ND1 to ND4 


observation, as it is not sensitive enough to measure sky 
radiances. The CCD version is about 1 0 4 times more sensi- 
tive (partly from the detector and partly from the addition of 
a lens to focus light on the optical fiber) and detects sky 
radiances easily. The CMOS head sensor is similar to the 
CCD version shown in Figure 2, but without the lens and the 
first filterwheel. Baffles are used to produce a purely geo- 
metrical field of view (FOV) of about 1.6° full width half 
maximum (FWHM), which is independent of wavelength. 
It uses an Avantes spectrometer with a 1024 pixel detector 
(Table 1). The Sun-and-Sky Pandora spectrometer system 
(Figures 1 and 2, Table 1) uses the head sensor as shown 
in Figure 2 and an Avantes spectrometer with a UV sen- 
sitive back-thinned 2048 x 16 pixel CCD. The neutral 
density filters (ND1, ND2, ND3, and ND4) in the first 
filterwheel give a 10,000:1 dynamic range, which com- 
bined with the exposure range 4-4000 ms gives an overall 
dynamic range of 10 7 needed for sun and sky observations 
at all latitudes and seasons. The 1.5° FOV of the Sun-and- 
Sky Pandora is produced by the lens and baffles with a 
small measured wavelength dependence. The Brewer FOV 
(2.6°) is larger than the Pandora FOV (~1.5°). Wavelength 
calibration and slit functions for the Pandora instruments 


Pandora Head Sensor 



Figure 2. Schematic of internal optics of the Pandora head sensor. 
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PANDORA DIRECT SUN SPECTRA 



Figure 3. The direct sun spectrum measured by Pandora 
using open hole (green), a U340 filter (blue), and the 
BP300 filter (red). 


are determined from lamp emission lines (Hg, Cd, Cu, In, 
Mg, Zn). Wavelength stability is validated during field use 
by an analysis of the solar Fraunhofer line structures. 

[ 9 ] Accurate measurements of TCO amounts in the 310 to 
330 nm range require a correction for the stray light 
effect, which is reduced by using UV band pass filters. 
Pandora instruments use two UV band pass filters, BP300 
(280-320 nm) and U340 (280-380 nm) with cutoffs at 
320 and 380 nm, respectively. The effect of the filters on 
the spectra is shown in Figure 3 compared to the spec- 
trometer response to the solar spectrum using an open 
hole in the filterwheel. A stray light correction (see section 3) 
is obtained using the signal on the detector from wavelength 
positions less than 290 nm where the incident sunlight is 
nearly zero. 

2.2. Pandora Total Ozone Spectral Fitting Algorithm 

[ 10 ] Here we analyze TCO amounts from Pandora oper- 
ated in direct-sun viewing mode. The spectral fitting algo- 
rithm uses laboratory measured absorption cross sections for 
each atmospheric absorber (i.e., 0 3 , N0 2 , S0 2 , HCHO, BrO 
for the 310-330 nm spectral region used here), a 4th order 
polynomial in wavelength to remove aerosols and Rayleigh 
scattering effects, and wavelength shift and squeeze func- 
tions. The shift and squeeze functions provide the best match 
to the measured spectrum compared to a solar reference 
spectrum containing the solar Fraunhofer line structure 
[ Herman et al., 2009]. High resolution O3 cross sections 
[Daumont et al., 1992], at a nominal effective ozone tem- 
perature of T = 225°K, were used in our retrievals. The 
effects of local and seasonal temperature changes on the O3 
cross sections and TCO retrievals are discussed in section 3. 

[ 11 ] The approach of using the independently measured 
Kurucz [2005] ground-based extraterrestrial spectrum, nor- 
malized to Atlas Susim-3, avoids the errors associated with 
performing a Langley extrapolation [e.g., Grobner and Kerr, 
2001; Bais, 1997] of the Pandora measured solar spectrum to 
zero air mass. The Langley extrapolation method, used by 
the Brewer, usually has to be made at a very clean atmo- 
spheric site such as Mauna Loa, Hawaii under the 


assumption that ozone is constant as a function of solar 
zenith angle (air mass). The resulting calibration can be 
transferred to other Brewer instruments, which then measure 
ozone relative to the transfer instrument. 

2.3. Brewer Total Ozone Retrieval 

[ 12 ] In its various versions, the Brewer spectrometer has 
been available since the 1980s for measurements of solar 
irradiance and zenith sky radiances, which can then be 
used to derive total column 0 3 , S0 2 , and N0 2 amounts [Kerr 
et al., 1981, 1985, 1988; Fioletov et al., 1998; Cede et al., 
2006b]. Most Brewers are single monochromators, but a 
number of systems are double monochromators with the 
advantage of a very low internal stray light fraction (<10~ 7 ) 
at the short UV-B wavelengths enabling high accuracy 
measurements even at high solar zenith angles (SZA). 

[ 13 ] The very low instrumental stray light makes the Mark 
III Brewer ideal for measuring Rayleigh scattered sky 
radiances in the UV spectral region except for the problem 
of polarization sensitivity related to the grating and to the flat 
quartz entrance window (Fresnel effect) covering the 
entrance port [Cede et al., 2006a], To eliminate both sources 
of polarization sensitivity so that accurate sky radiances can 
be obtained in the presence of aerosols, we have modified 
our Mark III Brewer instrument (#171) by introducing a 
depolarizer in front of the grating and by installing a curved 
quartz window such that the incident radiances observed by 
the spectrometer are always perpendicular to the window 
surface. Sky radiance measurements from Brewer #171 have 
been applied successfully to retrievals of 0 3 vertical distri- 
bution [Tzortziou et al, 2008]. For total column amount 
retrievals of ozone and other trace gases such as S0 2 and 
N0 2 , the modified Brewer system at Goddard has been used 
in direct-sun viewing mode. 

[ 14 ] The Brewer TCO retrievals are based on the Bass and 
Paur ozone absorption coefficients at a constant temperature 
of 228.3 K [Bass and Paur, 1985] that were modified by 
Kerr [2002]. TCO amounts were estimated here by using 
both the operational 4- and the 6-wavelength Brewer algo- 
rithms. In the operational algorithm the intensity of direct 
sunlight at the Brewer 4 longest wavelengths, 310.1, 313.5, 
3 16.8 and 320 nm, are used to calculate column ozone. A set 
of fixed weighting coefficients is applied, chosen to mini- 
mize the effect of S0 2 , temperature sensitivity, small 
wavelength shifts, and to suppress variations that change 
linearly with wavelength [Kerr et al., 1981, Kerr, 2002; 
Fioletov et al., 2005; Petropavlovskikh et al., 2011]. In the 
Brewer 6-wavelength algorithm, two additional wavelengths 
at 303 and 306 nm are used to calculate TCO, which 
enhances the retrieval sensitivity to ozone changes 
[Redondas and Cede, 2006]. However, the decrease in the 
signal at 303 and 306 nm with increasing SZA requires use 
of SZA-dependent weights in the 6-wavelength algorithm, 
such that at high SZA, 303 and 306 nm get low weights and 
the 6-wavelength algorithm is approximately reduced to 
the operational algorithm. Because of the relatively strong 
stray-light contribution in the 303-306 spectral region, the 
6-wavelength TCO retrieval is not recommended for mea- 
surements using a Brewer single spectrometer. 

[ 15 ] In the Brewer, there are 6 wavelengths in slit-mask 
mode, which are measured nearly simultaneously (0.1147 s 
per wavelength) reducing the influence of changing sky 
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conditions compared to classical wavelength scanning instru- 
ments (rotating grating). The Pandora system measures all 
wavelengths simultaneously, so that it has no “noise” intro- 
duced from changing sky conditions during a measurement. 
Although the Pandora instrument noise at a single wavelength 
is comparable to, or higher, than the Brewer, the noise effect in 
Pandora is minimized by using many more wavelength sam- 
ples (80 for the CMOS, 2x oversampled, and 160 for the CCD, 
4 x oversampled, spectrometers in the 3 10 and 330 nm spectral 
region), compared to 4-6 for the Brewer with no oversampling, 
and by time averaging (up to 200 individual measurements for 
the CMOS and 4000 for the CCD per 20 s). For example, at 
SZA = 70°, the Brewer noise at 3 10 nm (~0.7%) is better than 
the Pandora noise at 310 nm (~3%), while at 320 nm both 
instruments have about the same noise level (~0.4%). Simi- 
larly to the Brewer 6-wavelength algorithm, the Pandora 
algorithm weights the importance of the wavelengths inversely 
with the measured noise variance, so that as the SZA increases 
the importance of the short wavelengths decreases. Pandora 
stray light is comparable to single Brewer at about 1 x 10~ 5 
when normalized to 1 at the peak signal at each wavelength. 
However, the use of UV band-pass filters significantly reduces 
the stray light effect in Pandora (see section 3.1). 

2.4. Study Sites and Instrumentation 

[16] Greenbelt, Maryland: Measurements of TCO from 
our modified Mark III Brewer double monochromator (#171) 
and a collocated Pandora CMOS spectrometer (#9) were 
performed on the roof of a building, about 88 m above sea 
level, at GSFC in Greenbelt, MD (Lat = 38.993°N, 
Long = 76.840° W). The Brewer instrument was operating at 
GSFC from November 2000 to September 2010, providing a 
long time series of trace gas and aerosol amounts. Measure- 
ments from the Pandora CMOS instrument used in this study 
have been performed at this site since spring 2009. GSFC is a 
moderately polluted site, in the Baltimore-Washington DC 
metropolitan area. High ozone in the boundary layer is 
typically observed during April-September, particularly 
during the summer months [. Maryland Department of the 
Environment (MDE), 2008]. The site is located close to two 
major highway systems that are strong sources of NO x 
emissions resulting in high inhomogeneity and rapid changes 
in local tropospheric composition [Herman et al, 2009], 

[ 17 ] Cabauw, Netherlands: Total column ozone and trace 
gas retrievals by the Pandora #9 CMOS spectrometer were 
also performed at KNMFs Cabauw Experimental Site for 
Atmospheric Research (CAESAR) in the Netherlands 
(Lat = 5 1.970° N, Long = 4.926°E), during June-July 2009. 
Cabauw is a rural site, approximately 50 km from the North 
Sea and surrounded by densely populated and industrialized 
areas (33 km NE of the city of Rotterdam and 30 km SW of 
Utrecht). The climate in the area is characterized as a 
“moderate maritime” with a prevailing southwesterly circu- 
lation and precipitation spread over the seasons. Atmo- 
spheric conditions are characterized by large variability, 
ranging from very clean maritime to highly polluted air 
masses. Our measurements were part of the Cabauw Inter- 
comparison Campaign of Nitrogen Dioxide measuring 
Instruments (C1NDI campaign) [ Piters et al., 2012], 
Weather conditions during the campaign were mixed, with 
frequent changes in cloud cover, some rainy periods, and 
some early morning mist [e.g., Roscoe et al., 2010; Zieger 


et al, 2011]. Apart from NO 2 , other measured and inter- 
compared parameters included aerosol, HCHO, 0 2 0 2 , 
CHOCHO, and BrO. The campaign was held under the aus- 
pices of the European Space Agency (ESA), the international 
Network for the Detection of Atmospheric Composition 
Change (NDACC), and the EU Framework 6’s ACCENT- 
AT2 Network of Excellence and GEOMON Integrated 
Project. 

[is] Fairbanks, Alaska: Trace gas measurements by two 
Pandora instruments (the #9 CMOS and the #2 CCD sys- 
tems) and the Brewer double monochromator (#171) were 
performed in Fairbanks Alaska (Lat = 64.820°N, 

Long = 147.870°W) during March and April 201 1, as part of 
a large field campaign focusing on measurements of total 
column BrO, N0 2 and O3, and ozone vertical distribution at 
high latitudes. As part of this campaign, O3 and NO2 
amounts were obtained and compared to the Brewer double 
monochromator (#171). All instruments were located on the 
roof of a building at the University of Alaska Fairbanks 
campus, with a clear view of the morning sun horizon. The 
conditions during the campaign contained days that were 
mostly free of clouds (24, 25, 31 March and 1, 9, 1 1 April), 
some with light cloud cover (26, 27 March and 2 April), and 
the rest with moderate to heavy cloud cover. Both NCL and 
O3 data can be obtained by Pandora on days with light cloud 
cover, but with increased noise caused by the reduced 
number of measurements per 20 s. 

[ 19 ] Helsinki, Finland : Measurements from another Pan- 
dora CCD system (#21) are reported for Helsinki, Finland 
(Lat = 60.204°N, Long = 24.96 1°E). Observations have 
been performed as part of a campaign that started in Sep- 
tember 2011 to demonstrate the instrument’s performance at 
high SZA and to measure trace gases at high latitudes (NO 2 , 
O3, HCHO and BrO). The location of the instrument is on 
top of the Finnish Meteorological Institute (FMI), in an area 
of moderate automobile traffic with frequent maritime and 
occasional polluted air masses. 

3. Results and Discussion 

3.1. Ground-Based Instrument Intercomparison: 
Pandora and Brewer 

[ 20 ] To examine the Pandora performance over a relatively 
long data record, TCO retrievals from Pandora instrument #9 
were compared with measurements derived from the double- 
grating Brewer monochromator #171 at the GSFC site during 
June 2009 to September 2010 (Figures 4-6). Because of the 
higher frequency in Pandora measurements, the Pandora data 
were interpolated to Brewer measurement times, using only 
observations that were closer than 5 min apart. Pandora 
measurements were filtered for normalized root-mean square 
(RMS) of weighted spectral fitting residuals less than 0.05 
and estimated error in TCO less than 2 DU. Similarly, the 
Brewer measurements were filtered for standard error in TCO 
less than 2 DU. In addition, Brewer measurements were fil- 
tered for cloud and aerosols optical depth at 320 nm less 
than 1.0. 

[ 21 ] Comparisons at GSFC are shown between Pandora 
and Brewer for the 6-wavelength (6w) Brewer TCO retrievals 
(Figure 4a). Similar results were obtained for the operational 
(Op) 4-wavelength algorithm, although the Brewer and 
Pandora agreement was slightly better for the 6w-algorithm. 
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Figure 4. Comparison between TCO measurements from 
(a) Pandora (instrument #9) and the Brewer double-grating 
spectrometer (#171), collocated at the GSFC site, using the 
6-wavelength Brewer TCO algorithm, (b) Pandora (instru- 
ment #2) and the Brewer double-grating spectrometer 
(#171) in Fairbanks Alaska during March-April 2011. 

In both cases, the Pandora versus Brewer estimated least 
squares straight line fit was almost a perfect 1 : 1 within the 
standard deviation of the least squares fit (Slope = 0.987, 
Offset = 2.65 DU, R = 0.99, N = 2466, P < 0.0001 for the Op- 
algorithm; Slope = 0.991, Offset = 2.07 DU, R = 0.99, 
N = 2466, P < 0.0001 for the 6w-algorithm). The good 
agreement in accuracy over a wide range of ozone values 
(240-490 DU) suggests that the Brewer transferred Mauna 
Loa calibration and the absolute calibration of Pandora using 
the Kurucz and Atlas-3 SUSIM extraterrestrial spectrum are 
valid methods for measuring TCO. 


[ 22 ] A small bias of approximately 2 DU (or, 0.6%) was 
found between the Brewer and Pandora TCO measurements 
that was constant with slant column (SC) ozone amount over 
the full range of solar zenith angles (15-80°) observed dur- 
ing our measurements (Figure 5). For the Pandora compar- 
isons with the Brewer 6w-algorithm, scattering around this 
mean value of the bias was almost a constant ±10 DU (or 
±3%), independent of the SC amount or solar zenith angle 
(Figure 5a). Comparisons with the Brewer Op-algorithm 
(Figure 5b) resulted in somewhat larger scattering around the 
mean value at small SC ozone amounts, as would be 
expected, since the additional use of the 2 shortest wave- 
lengths (303 and 306 nm) in the 6w-retrieval improves the 
Brewer TCO estimates only at low SZAs or small SC (at 
high SZAs the 6w-algorithm is essentially reduced to the 
operational one). 

[ 23 ] Solar zenith angle dependence resulted in 0.4% 
change in TCO differences over the observed SZA-range 
(15°-80°) for the comparison between Pandora and the 
Brewer 6w-algorithm and 1% change in TCO for the com- 
parison with the Brewer Op-algorithm. Having the TCO 
differences almost independent of SZA (Figure 6) indicates 
that the stray-light method we applied to the Pandora mea- 
surements in the 310-330 nm region provided an adequate 
correction. However, the Brewer-Pandora differences 
showed some small seasonal dependence, with Pandora 
typically overestimating TCO during the winter months 
(square symbols in Figure 6) and underestimating TCO 
during the summer (circles in Figure 6). Average values of 


(A) 
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Brewer 6w SCO (DU) 


(B) 



Brewer Op SCO (DU) 

Figure 5. TCO residuals between Brewer and Pandora 
(instrument #9) as a function of Brewer Slant Column ozone 
(SCO), using (a) the 6-wavelength Brewer TCO algorithm 
and (b) the operational Brewer TCO algorithm. The practically 
zero slant column 0 3 dependence suggests a negligible remain- 
ing stray light effect after the applied simple SL correction. 
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SZA (°) 


Figure 6. TCO residuals between Brewer (6w-algorithm) and Pandora (instrument #9) as a function of 
solar zenith angle (SZA). Results are shown in different symbols for different seasons (triangles for spring, 
MAM; circles for summer, JJA; diamonds for fall, SON; squares for winter, DJF). The liner regression on 
the data is shown as solid line and reflects the almost negligible SZA-dependence of the Brewer-Pandora 
residuals. 


Brewer-Pandora differences were +3.6 DU (s.d. = 3.8 DU) 
or +1.15% during summer, and —4.74 DU (s.d. = 4.3 DU) or 
— 1.32% during winter. Better agreement was observed 
during the spring and fall months with average TCO differ- 
ences between the two instruments of 1.3 DU (s.d. = 4.2 DU) 
and 2.4 DU (s.d. = 3.6 DU), respectively. This small seasonal 
variability could be associated with the temperature depen- 
dence of the 0 3 cross sections, assumed constant (T = 225°K) 
in our standard retrieval. Our sensitivity studies show that the 
retrieval of TCO has a temperature dependence of 1% per 
3°K, with an underestimation in temperature (e.g., during 
summer) resulting in an underestimation of TCO values (i.e., 
positive Brewer-Pandora differences). Application of a tem- 
perature correction by approx. ±5°K, consistent with the 
TOMS version 8 temperature climatology for the specific 
latitude and TCO values [McPeters et al., 2007], consider- 
ably reduces this small (<1.5%) seasonal dependence in TCO 
differences between the two instruments. 

[ 24 ] Almost no SZA dependence was found for the com- 
parison between the Brewer and Pandora instruments at the 
high-latitude site in Fairbanks, Alaska (Figures 4b, 8d) over 
the 58-80° SZA range observed by Brewer. TCO measure- 
ments were almost identical between the two Pandora 
instruments, the CMOS (#2) and CCD (#9) spectrometer 
systems (Figure 8d). Linear regression between Brewer and 
Pandora #2 resulted in a slope of 0.979 and strong correla- 
tion (R = 0.99). Flowever, an offset of 16.55 DU was found 
between the two instruments, with Pandora systematically 
lower by 2^1% relative to the Brewer (Figures 4b and 8d). 
This is larger than the 2 DU offset found during the much 
longer GSFC Brewer-Pandora time series. 

3.2. Comparison Between Pandora and Aura OMI 

[ 25 ] Close agreement was found between Aura OMI and 
our ground-based measurements of TCO over all four study 
sites (Figures 7, 8, 10, 11). For comparisons with satellite 
observations, Brewer and Pandora data were averaged over a 
±l-h window of the OMI overpass time. Ground-based 
measurements were filtered to only include data with 
uncertainty in TCO <2 DU and normalized RMS of 


weighted spectral fitting residuals <0.05. OMI data were 
filtered to keep only those data where the distance from 
GSFC was <50 km and the OMI radiative cloud-fraction 
was <0.2. 

[ 26 ] Both Brewer and Pandora data agreed closely with 
the OMI near-overpass at GSFC, with an estimated slope not 
significantly different from 1 (1.001 and 1.009 respectively, 
Figures 7a and 7b). Both Pandora and Brewer slightly 
overestimated TCO relative to OMI by 0.7% and 0.3% 
respectively, with residuals in the range ±7%. Satellite and 
ground-based (Pandora or Brewer) TCO differences showed 
some dependence on OMI cross-track position (CTP) and 
the distance between the station and the OMI CTP, with the 
largest TCO differences (larger than 10%) associated in 
almost all cases and all sites with extreme CTP and distance 
from station larger than 25 km. 

[ 27 ] Very good agreement between satellite and Pandora 
retrievals was also observed at Helsinki, Cabauw and Fair- 
banks (Figures 7c-7e), although the Pandora-OMI match- 
ups were limited at these locations due to the shorter periods 
of ground-based instrument deployment. In all cases, linear 
regression between OMI and Pandora retrievals resulted in a 
slope close to 1 (0.95-1.02) and strong correlation coeffi- 
cients (R within 0.97-0.99). Residuals between ground- 
based and satellite data at Cabauw ranged between ±2%, 
with an average difference between Pandora and OMI of just 
0.3 DU (s.d. = 3 DU). Pandora underestimated TCO relative 
to OMI at Fairbanks, Alaska, on average by 1% (—4 DU 
average Pandora-OMI difference, s.d. = 5 DU), with TCO 
residuals between —12 and 3 DU. Differences between OMI 
and Pandora retrievals were within ±5% at the FMI station 
in Helsinki, Finland, with Pandora underestimating TCO, on 
average by 0.7% (Figures 7c-7e). 

3.3. Temporal Variability in TCO From Pandora, 
Brewer, and OMI 

[ 28 ] The frequent and continuous measurements by Pan- 
dora from sunrise to sunset allowed us to study temporal 
dynamics in ozone at our study sites (day-to-day and diurnal 
variation), and revealed short-term changes in TCO amounts 
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Figure 7. Comparison in TCO between AURA-OMI and (a) Brewer 6w-algorithm at the GSFC site 
(N = 90); (b) Pandora #9 at the GSFC site (N = 99); (c) Pandora #21 at Flelsinki (N = 32); (d) Pandora 
#9 at Cabauw (N = 12); and (e) Pandora #9 at Fairbanks (N = 18) (similar results for Pandora #2 at Fair- 
banks). In all cases, P <0.0001. 


not possible to capture by sun synchronous satellite obser- 
vations, like Aura OMI, alone. 

[ 29 ] Figure 8a shows comparisons between Pandora, 
Brewer and OMI over GSFC, for four example cases cov- 
ering different seasons during the 2009-2010 period when 
Brewer and Pandora were collocated at GSFC. Similar 
comparisons between Pandora and OMI (and Brewer where 


available), are shown for Cabauw, Helsinki and Fairbanks in 
Figures 8b-8d, for four example cases covering the period of 
instrument deployment at these sites. In general, the diurnal 
TCO patterns depicted by our two ground-based instruments 
Pandora and Brewer were in very close agreement, both at 
GSFC and in Fairbanks, with few exceptions. On 4 May 
2010, the morning data (9:00 to 11:00 UT) at GSFC show 
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8 . A comparison of total column ozone amounts between Pandora instruments (diamonds), 
6 w-algorithm (circles), and OMI (squares) at (a) GSFC, (b) Cabauw, (c) Helsinki, and (d) Fair- 


very good agreement. After that, the Brewer data become 
sparse, while the Pandora data remain at high precision. The 
usual reason for this type of behavior in the Brewer is that 
the aerosol variation was significant over the period (0.7 s) 
of the Brewer sequential wavelength measurement, causing 
significant O 3 variation between successive retrievals used 
in forming an average value (and standard error in estimated 
Brewer TCO > 2 DU). Aerosol variation has little effect on 
the Pandora TCO retrieval, plus Pandora can average more 
measurements together to minimize the noise in a single data 
point. Consistent with our results on the SZA-dependence of 
Brewer-Pandora TCO differences (discussed in section 3.1), 
the observed overall good agreement between the two 
instruments throughout the day suggests minimum residual 
stray light effects after stray light correction. This shows that 


Pandora can be used for satellite validation and ozone 
characterization under a wider variety of atmospheric con- 
ditions compared to the Brewer spectrometer, while being of 
comparable accuracy under clear-sky conditions. 

[ 30 ] Variation in satellite retrieved total column ozone 
amounts over the approximately 1.5 h between two succes- 
sive OMI orbits ranged from less than 1 DU (i.e., 24 June 
2009, at Cabauw) to more than 19 DU (i.e., 4 May 2010, at 
Goddard) for the days shown in Figure 8. The large increase 
in TCO observed by OMI over GSFC on May 4th 2010 was 
also observed by Pandora and Brewer, indicating a real 
change in ozone amount. Such sharp (10-20 DU) spikes and 
dips in TCO have been observed in other days by both OMI 
successive overpasses and the higher resolution Pandora, 
and are mostly associated with changes in stratospheric 
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Figure 8. (continued) 


ozone content due to weather systems (e.g., pressure chan- 
ges, passage of a cold front with high ozone content, or 
intrusion of low-ozone air from lower latitudes). Aura-OMI 
observations of the ozone field over the Eastern U.S. on May 
4th 2010 revealed that the sharp temporal change in TCO we 
observed at Goddard that day was indeed due to the passage 
of a front transporting air masses with higher TCO from the 
north to the south and over the GSFC area (Figure 9). 

[ 31 ] Differences in retrieved TCO between two successive 
OMI overpasses, however, do not always reflect real chan- 
ges in ozone content over a specific location, as this could 
also be the result of differences in the satellite FOV, large 
changes in distance between the station and the satellite 
cross-track position CTP (i.e., by more than 50 km), or rapid 
changes in cloud fraction and aerosols between the two 
overpasses. For example, on 27 April 2010 (data not shown 
here) the change in OMI retrieved TCO between two 


successive satellite overpasses at GSFC was as high as 
47 DU, with a change in station-CTP distance, however, of 
more than 80 km. Both Brewer and Pandora at GSFC mea- 
sured a change in TCO of less than 10 DU between the two 
OMI overpass times, suggesting a relatively large influence 
of horizontal atmospheric heterogeneity on the OMI TCO 
retrieval. 

[ 32 ] For the cases shown in Figure 8, diurnal TCO vari- 
ability as high as 40-50 DU was observed by our Pandora 
instruments (e.g., 30 April 2010 at GSFC, 2 April 2011 at 
Fairbanks). Such temporal changes in ozone occur occa- 
sionally at these mid- to high-latitude sites. An average 
amplitude of diurnal variation in TCO of 14 DU has been 
measured by Pandora over Goddard during the period May 
2009-April 2011, while temporal changes as strong as 40 to 
50 DU have been recorded on certain days by both Brewer 
and Pandora at this site. A diurnal variability of similar 
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Figure 9. Aura-OMI TCO over the Eastern U.S. for two successive orbits on May 4, 2010 (orbits L2- 
OMTO3_20 1 0m0504t 1 627-o30864_v003- 20 1 2m033 1 122 1 500.he5 and_L2-OMT03_2010m0504tl806- 
o30865_v003-2012m0331t221553.he5). The OMI images reveal the passage of a front moving from North 
to South and transporting air masses with higher TCO over the GSFC area. This spatial heterogeneity in 
TCO was reflected in the TCO time series measured during this day by both Brewer and Pandora at the 
GSFC site (see Figure 8a). Aura-OMI satellite images provided by Gordon Labow. 


magnitude was observed in both Cabauw and Fairbanks, 
with an average change in TCO during the day of 16 DU and 
20 DU, respectively. Diurnal variation with amplitude as 
strong as 25 to 30 DU was measured in several occasions 
over the 2 months (June-July) of the CIND1 campaign. TCO 
changes as large as 30 to 46 DU were observed in Fairbanks, 
by both Brewer and Pandora during March 2011. Observed 
variability in TCO was lower at Helsinki, Finland, where 
average TCO change during the day was 12 DU during 
September to November 2011, with a maximum observed 
amplitude of diurnal variation of 33 DU. Observations at our 
high-latitude site in Helsinki during days when TCO 
remained almost constant throughout the day, revealed that 



Time (UT) 


Figure 10. Comparison between Pandora (open diamonds) 
and OMI (filled square) for 24 October 2011, a typical day 
when TCO remained constant throughout the day, at the 
high latitude site Helsinki, in Finland (Lat = 60.2037°, 
Long = 24.9612°). The difference between Pandora and 
OMI TCO was 7.26 DU, or 2.8%. Pandora measurements 
are shown during the day (time on x axis and SZA on right 
y axis) The effect of residual stray light correction in Pan- 
dora becomes apparent only at SZAs larger than 80 degrees. 


the effect of stray light correction on the Pandora retrievals 
becomes apparent only at SZAs larger than 80 degrees 
(Figure 10). In terms of the slant optical path that includes 
total ozone dependent extinction, these results saw that 
Pandora TCO retrievals are accurate up to slant column 
ozone (SCO) amounts of 1,500 DU (=260 DU * AMF(80°), 
where AMF is the air mass factor). This is consistent with 
results at GSFC showing good agreement between Brewer 
and Pandora retrievals up to slant column ozone amounts of 
approximately 1,500 DU (Figure 5). Thus, in addition to the 
tight filtering of the data that removed all outliers (normal- 
ized root-mean square of weighted spectral fitting residuals 
less than 0.05 and estimated error in TCO less than 1 DU), 
measurements at SCO larger than 1500 DU were excluded 
for the above estimates of diurnal variability at our four 
study sites. 

[ 33 ] The long time series of TCO measurements at 
Goddard was used to estimate seasonal variability in ozone at 
this midlatitude site (Figure 1 1). Estimates of monthly means 



Month 

Figure 11. Mean monthly TCO values at GSFC measured 
by Brewer (circles), Pandora instrument #9 (diamonds) and 
OMI (squares). The TCO zonal-mean monthly mean clima- 
tology from OMI/MLS at 39° latitude [Ziemke et al, 2011] 
is also shown for reference (gray line). 
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showed considerably higher ozone amounts and high vari- 
ability during late winter and early spring, when monthly 
mean TCO was between 320 and 380 DU. A second smaller 
maximum was observed during July (~320 DU), possibly 
due to the particularly high tropospheric ozone amounts 
during summer in this area [MDE, 2008; Ziemke etal, 2011]. 
Total column ozone dropped during the fall months, with 
monthly mean TCO in the range 275-300 DU. These results 
are consistent with the TCO zonal-mean monthly mean cli- 
matology derived from OMI and MLS (Microwave Limb 
Sounder) ozone measurements (2004-2010) at 39° latitude 
[. Ziemke et al., 2011]. Comparisons of TCO monthly avera- 
ges between Pandora and OMI showed differences <6%, 
while Brewer and OMI showed differences <8%. Although 
these are observations at individual sites and over relative 
short time periods, Pandora measurements at our mid- to 
high- latitude sites (Cabauw at 54°N, Helsinki at 60°N and 
Fairbanks at 65°N) showed seasonal and latitudinal patterns 
in TCO that were generally in agreement with climatology. 
That is, an increase in TCO would be expected, based on 
climatology, with increasing latitude in the Northern Hemi- 
sphere, with higher values in early spring and lower values in 
fall [Ziemke et al, 2011], Total column ozone amounts were 
particularly high at our highest latitude site in Fairbanks 
Alaska (Figure 8). Monthly mean TCO was 426 DU (s. 
d. = 19 DU) in March and 420 DU (s.d. = 22 DU) in April 
2011, approximately 10% higher than the 60-65° latitude 
monthly mean climatology from OMI and MLS. Total col- 
umn ozone at the midlatitude Cabauw site during the CINDI 
campaign was on average 328 DU (s.d. = 15 DU) in June and 
315 DU (s.d. = 17 DU) in July, very similar to the summer 
monthly TCO means at the midlatitude GSFC site. Lower 
ozone values were observed during fall 2011 in Helsinki, 
with monthly TCO of 284 DU (s.d. = 19 DU), 277 DU (s. 
d. = 31 DU), and 257 DU (s.d. = 14 DU), in September, 
October and November, respectively, consistent with TCO 
climatology in the 55-65° latitude zone. 

4. Summary 

[ 34 ] The Pandora and Brewer TCO measurements at the 
GSFC site were found to be in close agreement (0.6% bias) 
throughout the day, including near sunrise and sunset con- 
ditions with solar zenith angles up to 80°, or SCO amounts 
up to 1,500 DU. The good agreement in accuracy over a 
wide range of ozone values (240-490 DU) suggests that the 
Brewer transferred Mauna Loa calibration and the absolute 
calibration of Pandora using a combination of the indepen- 
dently measured Atlas-3 SUSIM and Kurucz extraterrestrial 
spectra are in good agreement. The practically negligible 
SZA dependence of the TCO residuals between Pandora and 
the well-calibrated double grating Brewer spectrometer #171 
that was observed at both GSFC and Fairbanks and for both 
Pandora CMOS (#2) and CCD (#9) systems, shows that the 
stray-light correction method applied to Pandora provided an 
adequate correction. This successful ground-based instru- 
ment intercomparison suggests that the Pandora direct-sun 
spectral fitting technique can be reliably used for satellite 
validation (e.g., Aura-OMI and NPP-OMPS), and can be 
applied to characterization of total column ozone dynamics 
under a wider variety of atmospheric conditions compared to 
direct sun observations from the Brewer spectrometer, while 


being comparable under clear-sky conditions. Zenith sky 
observations from both instruments could further extend 
TCO retrievals under cloudy conditions [. Fioletov et al, 
2011]. Pandora and Aura-OMI TCO retrievals were found 
to be in very good agreement across all sites, GSFC, 
Cabauw, Helsinki and Fairbanks, when the OMI cross track 
position was at a distance of less than 50 km from the Pan- 
dora location, and OMI-measured cloud fraction was less 
than 0.2. In all cases, linear regression between OMI and 
Pandora retrievals resulted in a slope close to 1 (0.95-1 .02) 
and strong correlation coefficients (R within 0.97-0.99). 
Residuals between ground-based and satellite data ranged 
between ±7%, with Pandora slightly overestimating TCO at 
GSFC (by 0.7%), while slightly underestimating TCO at 
Fairbanks (1%) and Helsinki (0.7%). 

[ 35 ] Our Pandora ground-based measurements at four sites 
extending from 39°N to 65°N latitude, revealed seasonal and 
latitudinal patterns in TCO that were generally in agreement 
with OMI-MLS climatology. Overall, we observed increas- 
ing TCO with increasing latitude, with higher values in early 
spring and lower values in fall. A second smaller maximum 
was observed during July at GSFC, most likely associated 
with the particularly high tropospheric ozone amounts typi- 
cally observed in this area during summer. Considerable 
diurnal variability, on some days as high as 40 to 50 DU, 
was capture by the frequent and continuous Pandora mea- 
surements at these mid- to high- latitude sites. Short-term 
(hourly) temporal changes in ozone amounts cannot always 
be captured by OMI or other satellite instruments in sun- 
synchronous orbits with 1 to 3 daily overpasses. Such a 
variation in TCO, of the order of 10 to 15%, is especially 
important if it involves an ozone reduction near noon hours, 
since it causes a 10 to 15% increase in erythemal UV irra- 
diance, which is an important factor in human health. Such 
variations are not usually included in the UV index forecast. 
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